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This paper presents a comprehensive study of the fundamental quasinormal modes of all Stan- 
dard Model fields propagating on a brane embedded in a higher-dimensional rotating black hole 
spacetime. The equations of motion for fields with spin s = 0,1/2 and 1 propagating in the 
induced-on-the-brane background are solved numerically, and the dependence of their QN spectra 
on the black hole angular momentum and dimensionality of spacetime is investigated. It is found 
that the brane-localised field perturbations are longer-lived when the higher-dimensional black hole 
rotates faster, while an increase in the number of transverse-to-the-brane dimensions reduces their 
lifetime. Finally, the quality factor Q, that determines the best oscillator among the different field 
perturbations, is investigated and found to depend on properties of both the particular field stud- 
ied (spin, multipole numbers) and the gravitational background (dimensionality, black hole angular 
momentum parameter). 

PACS numbers: 04.30.Nk,04.50.+h 
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I. INTRODUCTION 

The introduction, during the recent years, of gravi- 
tational theories that postulate the existence of addi- 
tional spacelike dimensions in nature 0, Q] have radi- 
cally changed our notion of the universe and of the struc- 
ture of the spacetime in which we live. The concept of 
the brane - a (3+l)-dimensional hypcrsurface inspired 
by D-branes encountered in string theory - was intro- 
duced, that was embedded in a higher-dimensional space- 
time, the bulk. Ordinary Standard Model (SM) particles 
(scalars, fermions and gauge bosons) are restricted to live 
on the brane while gravitons can propagate both on and 
off the brane. In this set-up, the accurately observed 
properties of SM fields are only marginally altered, nev- 
ertheless new effects coming from the embedding of our 
4-dimcnsional world in a higher-dimensional spacetime 
could still be observed, despite the fact that the bulk 
itself remains strictly non-accessible to brane observers. 

One such effect could be the creation of microscopic 
black holes on ground-based accelerators during the colli- 
sion of highly energetic particles with center-of-mass en- 
ergy y/s > M* :3l], where M* denotes the fundamental 
Planck scale of the higher-dimensional gravitational the- 
ory. This scenario can be realized in the context of the 
theory with Large Extra Dimensions p| that allows for 
a fundamental gravity scale as low as 1 TeV. As a re- 
sult, experiments performed at next-generation particle 
colliders could easily witness effects coming from trans- 
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planckian collisions, such as the creation of microscopic, 
higher-dimensional black holes that are centered on our 
brane but extend off the brane as well. The creation of 
these black holes, and thus the existence of additional 
dimensions, could be confirmed through the detection of 
the emitted Hawking radiation which is expected to be 
their most prominent observable effect 

Nevertheless, the creation of these brane black holes 
might trigger additional observable effects, such as the 
detection of the spectrum of the so-called quasi-normal 
modes (QNMs) [H, ||| . These modes arise as the result 
of an external perturbation of a black hole background, 
either through the addition of a field or by perturbing 
the metric itself. When this happens, the gravitational 
system enters a phase of damping oscillations, with the 
frequency of the field consisting of a real part wr c , that 
drives the field oscillations, and of an imaginary part 
LUi m , that causes the simultaneous damping of these oscil- 
lations. Quasi- normal modes with small u>i m , and thus 
long damping time, can dominate the spectrum at very 
late times after the initial perturbation. 

As the detection of the spectrum of QNMs could be 
considered as direct evidence for the existence of black 
holes in our universe, a variety of black-hole backgrounds 
and auxiliary propagating fields have been studied over 
the years both in 4 dimensions [j| and higher dimensions 
0. Unfortunately, no experimental detection has so far 
confirmed these theoretical studies. This could be due 
either to the notorious elusiveness of gravitons or to the 
fact that simply no black hole exists, as yet, in the vicin- 
ity of our neighbourhood. If, however, the scenario of 
the creation of microscopic black holes can be realised at 
the next-generation colliders, a unique opportunity for 
detecting the QN spectrum is offered. Since we have not 
yet successfully detected any gravitons, the obvious can- 
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didates would be instead the different species of SM fields 
that are not only abundant on our brane but also much 
easier to detect. 

In we presented a comprehensive study of the QN 
spectra of all SM fields (scalars, fermions and gauge 
bosons) living on a 4-dimensional brane, that was em- 
bedded in a variety of higher-dimensional, spherically- 
symmetric black-hole backgrounds. The cases of brane- 
localised fields propagating in a Schwarzschild, Reissner- 
Nordstrom and Schwarzschild- (Anti) de Sitter induced 
gravitational background were studied in detail, and the 
effect of various parameters, such as the dimensionality 
of spacetime, black hole charge and bulk cosmological 
constant, on the QN spectra was examined. It was found 
that an increase in the number of transverse-to-the-brane 
spacelike dimensions resulted in the faster damping of 
all fields living on the brane, independently of the ex- 
act type of spherically-symmetric projected background. 
In the presence of a black hole charge, the QN spectra of 
all fields resembled more the ones of Z?-dimensional fields 
rather than the ones of 4-dimensional, while the presence 
of a bulk cosmological constant did not lead to significant 
modifications to the QN spectra compared to the ones of 
purely 4-dimensional fields. A recent paper ^(j follows 
the same line of research by studying the quasinormal 
modes of brane-localised black holes and their thermo- 
dynamic properties. 

In the present work, we extend our previous analy- 
sis by considering higher-dimensional black-hole back- 
grounds with non-vanishing angular momentum. During 
the high-energy collisions of elementary particles result- 
ing in the creation of black holes, it is unnatural to ex- 
pect that only head-on collisions, leading to spherically- 
symmetric black holes, would take place. Collisions with 
a non-vanishing impact parameter are much more likely 
to occur, and, in addition, it is for these collisions that the 
black-hole production cross-section is maximised 0, Q . 
Therefore, microscopic rotating black holes should be the 
most generic situation, and the effect of the angular mo- 
mentum of the black hole on the QN spectra of brane- 
localised fields needs to be investigated. An earlier paper 
[TT| that appeared in the literature studied the QN spec- 
tra of brane-localised SM fields only for the particular 
case of a 5-dimensional rotating black hole. Here, we 
extend this work by presenting a comprehensive analy- 
sis of the QN spectra of fields with spin s = 0, 1/2 and 
1 for arbitrary dimensionality of spacetime and angular 
momentum of the black hole. 

As in our previous analysis, in order to ignore quantum 
corrections, we will be assuming that the produced black 
hole has a mass Mbh that is at least a few orders of 
magnitude larger than the fundamental scale of gravity 
M*. Also the brane self-energy can be naturally (i.e. in 
the context of a theory that solves the hierarchy problem) 
assumed to be of the order of M*; therefore, it is much 
smaller than Mbh, and its effect on the gravitational 
background can also be ignored. 

In section II, we present the theoretical framework for 



our analysis and the equations of motion for the SM fields 
propagating in the brane background. Section III out- 
lines the numerical techniques employed for the integra- 
tion of both the radial and angular part of the equation of 
motion of a field propagating in an axially-symmetric in- 
duced background. Our numerical results are presented 
in Section IV, where the effect of the dimensionality of 
spacetime, angular momentum of the black hole and spin 
of the particle on the QN spectra is investigated. We fin- 
ish with your conclusions in Section V. 



II. MASTER EQUATION FOR PROPAGATION 
OF FIELDS ON THE BRANE 



The line-element describing a higher-dimensional, ro- 
tating, neutral black hole is given by the Myers-Perry 
solution As mentioned above, in this work we will 
concentrate on the propagation of SM fields, i.e. fields 
with spin s = 0, 1/2 and 1, on the induced-on-the-brane 
gravitational background. This background is given by 
the projection of the higher-dimensional one onto the 
brane by fixing the values of the additional angular coor- 
dinates that describe the (D — 4) extra spacelike dimen- 
sions 0,0. After the projection, the brane background 
assumes the form H 



ds 2 
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dr 
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2a/! sin 9 



a 2 /x sin 2 I 

VJ r D-5 



dt dip 



-dr 2 



(1) 
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where 



A = r 2 + a 2 



r D-5 



and S = r 2 + a 2 cos 2 9 . (2) 



The parameters /i and a are related to the mass and an- 
gular momentum, respectively, of the black hole through 
the definitions I12I 



(D - 2) tt^" 1 )/ 2 



J 



D 



■ M B n a . 



with k 2 d = 8ttG = 8tt/M® 2 the D-dimensional New- 
ton's constant. We should note here that the higher- 
dimensional black hole is assumed to have only one non- 
vanishing component of angular momentum, about an 
axis in the brane. This is due to the simplifying assump- 
tion that the particles that created the black hole were 
restricted to live on an infinitely- thin brane [lj, therefore, 
during collision they had a non- vanishing impact parame- 
ter only on a 2-dimensional plane along our brane. Also, 
one should observe that, as in the case of spherically- 
symmetric black holes, the induced-on-the-brane back- 
ground has an explicit dependence on the total number 
of dimensions including the ones transverse to the brane 
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For the study of the quasi-normal modes of all SM 
fields living on the brane, we need first to derive their 
equations of motion for propagation in the induced back- 
ground (JIJ. For this, the following factorised ansatz 

*.(t, r, 8, if) = e- 1 ^ A~ s P s (r) S™ e {6) , (4) 

is employed for a field with spin s, where S™ e (9) are 
the spin- weighted spheroidal harmonics ^5|. Then, the 
use of the Newman-Penrose |16L Il7j formalism may lead 
to a "master" equation describing the propagation of all 
species of SM fields (scalars, fermions and gauge bosons) 
on the brane. This equation has been derived in |illl8j. 
and, as in the case of a purely 4-dimensional Kerr back- 
ground 0, it leads to two decoupled equations for the 
radial 



A s 



dr 



A 1 " 



dP s 
dr 
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+ Us co r - A )P S = , (5) 



and angular part 
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dS 
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sm ( 
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2atos cos 9 + s 



s 2 cot 2 I 
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of the field, respectively. In the above, we have used the 
definitions 

K = (r 2 +a 2 ) to— am , A = \+2s+a 2 oj 2 —2amio . (7) 

As in 4 dimensions [l^, the above equations hold for 
the upper component s = \s\ of all SM fields with spin 
s = 0, 1/2 and 1, but in this case they describe their 
propagation not in a purely 4-dimcnsional background 
but on a brane embedded in a higher-dimensional Kerr 
background. 



where 



2 2,2 / 2 . 2 

p =r +a —am/uj — r +a 



(9) 



and the metric function A(r) is given in Eq. (J2J. Then, 
Eq. © takes the form 



\ dr 2 ' J V dr* 

where we have defined 



icj)Y s -QY s = 0, (10) 



P = s 



/4rA A' 



(11) 



and 



A - (2s- 1) 



A 

Ar 



A'r 

+ (2a-3) — - (s - 1) — 
P P 



By defining further, as in |20(, 



Y s = hZ s + 2iuj ( iuj ) Z s , 



(12) 



(13) 



Eq. i|l(J|) can now be written as a one-dimensional Schro- 
dinger equation 



if 
dr 



— +lu 2 )Z s = V s Z s , 



with the effective potential given by 



~ 2 A . 2 d (A 
X ~ a 7 TlaP Tr \7 



(14) 



(15) 



for spin-1 particles, 



A. Alternative form: One-Dimensional wave 
equation 

As in the case of the study of QNMs of SM fields living 
on a projected-on-the-brane spherically-symmetric black- 
hole background 9], we will find it convenient to rewrite 
the radial equation in the form of a one-dimensional 
Schrodinger (or, wave-like) equation. Here, we gener- 
alise both the analysis of [20] performed for the case of a 
purely 4-dimensional Kerr-like background, as well as the 
one of where a £>-dimensional spherically-symmetric 
background was projected on the brane. To this end, 
we define a new radial function and a new ("tortoise") 
coordinate according to 



.2(8-1/2) Y 



dr* 
dr 



£1 
A 



(8) 



A 

y s =i/2 = — 
p 



dr 



(16) 



for spin- 1/2 particles, and 



V s -_ 



P 



A + 



A 3Ar 2 



A'r 

~7 



(17) 



for spin-0 particles. The QN modes for gauge bosons, 
fermions and scalars living on the brane in the vicin- 
ity of a Kerr-like induced background can be found by 
solving Eq. i|Ti|) upon choosing the appropriate effective 
potential. Here, we will undertake this task and compute 
the QNMs of brane localised SM fields through numeri- 
cal analysis the details of it are explained in the next 
section and our results are presented in the following one. 
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III. NUMERICAL ANALYSIS 

In order to calculate the QNMs of brane-localized 
fields, we need to solve the radial equation or its 
equivalent form (|14|l . To this end, we need the expression 
of the angular eigenvalue A, that appears in the radial 
equation. For a vanishing black hole angular momen- 
tum, i.e. a = 0, the angular eigenvalue is simply given 
by A = £(£+1) - s(s+ 1), where £ = s, s + 1, s + 2 . . . are 
(half)integer numbers; therefore, the angular equation 
can be altogether ignored for the purpose of calculating 
QN modes. However, in the case of an axially-symmetric 
black-hole background, the angular eigenvalue cannot be 
written in a closed form, and its exact value can be found 
only numerically [2l| by using the angular equation JSJ). 

For a ^ 0, the spectrum of A for a given value of u> is 
a set of complex numbers that cannot be easily put into 
correspondence with a set of angular numbers I. In order 
to find the exact value of X(u>), and at the same time the 
solution to Eqs. J3IHJ), we follow a technique inspired by 
perturbation theory: we start with the known solution 
of the angular equation for a = (i.e. with the spin- 
weighted spherical harmonics), and increase a by a small 
amount; then, for a given energy w, we search for the 
solution of Eqs. (JSJ HJJ) that is regular and closest to the 
original one; this solution is characterized by a unique 
value \{oS). By repeating the same analysis for different 
values of (a,o>), the spectrum of A(w), as well as the 
solution to the radial and angular equations, can be found 
for a wide range of the angular momentum parameter a. 

The quasinormal modes follow from the numerical in- 
tegration of the radial equation J5J after imposing the 
following boundary conditions |3, |(j 



P, 



P, 



C+ exp(io;r)/r , 
C.(r-r h )- iK , 



as 



as 



oo . 



rh ■ 



where 



ur h (rl + a 2 ) - mar h 
(n-l)(rl + a2) + 2rl 



(18) 
(19) 



(20) 



and thus correspond to outgoing waves at spatial infinity 
and to ingoing waves at the event horizon r^. The black 
hole horizon is derived by solving the equation A(r) = 0. 
Unlike the 4D Kerr black hole for which there is an inner 
and outer solution for rh , for D > 5 there is only one real, 
positive solution to this equation, satisfying the relation 



„£>-3 



/x/(l + a 2 ), where we have defined a* = a/rv 



The radial equation JSJ) can be solved by using the 
continued fraction method p2l |. For this, the remaining 
regular singular points of Eq. (JSJ need to be identified. 
For D = 5, there is only one more singularity at r = — rv,. 
That is why the Frobenius series has a form similar to 
the one of the ordinary 4-dimensional Kerr background 



= exp^r) ( r_rh\- lK y ^ ( r_Th 
r + rh \r + r h J ^ \r + r h 



(21) 



For D > 5, there are D — 3 additional singularities, one of 
them being at r = 0. Fortunately, for D < 9, |1 — r^/r| > 
1 for all of them, thus the appropriate Frobenius series 
in this case has the following form 



Ps(r) 



exp(iwr) / rh \ ~ tK 



0-t) 2>0- t) ■ <22) 



i=0 



For D > 9, some of the singularities appear in the 
unit circle and one has to continue the Frobenius series 
through some midpoints. This technique was recently 
developed in [23]. 

Substituting the above series into Eq. (0), one can 
obtain a (2D — 5)-term relation for the coefficients hi 

min(2D—6,i) 

c£f- 5 V, A) &*_,=(), fori>0. (23) 

3=0 

Each of the functions c^f ^ (u>, A) has an analytical form 
in terms of s, i, u) and A, for a specific dimensionality D 
and each integer < j < 2D — 6. 

We now decrease by one the number of terms in the 
relation 



min(k.i) 

£ c^ +1 V, A) fa^ = , 

3=0 

i.e. find the c^(lu,X) that satisfy the equation 

min(k — 
3=0 

For i > k, we can rewrite the above as 

O+l)/ \N 



(24) 



(25) 



4-i,i-i( w > A ) j^i 

Subtracting from J23J we find the relation (|2*S)l ex- 
plicitly, thus we obtain: 

eg 5 (w, A) = cg +1) (w, A), for j — 0, or i < k, 



'£cf\ i _ 1 (uj,\)b i ^ = 0. (26) 



C ^( w ,A) = c^(c,A) 



otherwise. This procedure is called Gaussian elimina- 
tions, and allows us to determine numerically the coeffi- 
cients in the three-term relation 



c§ 6< + 45 «N-l + 45 &<-2 - 0, for i > 1 
4 1 fti + 4 i = 0, 



(27) 



for given lu (and therefore A) up to any finite i. The 
complexity of the procedure is linear with respect to i 
and k. 
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The requirement that the Frobenius series be conver- 
gent at spatial infinity implies that 



= c 



(3) 
1,1 



„(3) (3) (3) (3) 
-0,1 L 2,2 L 0,2 L 2,3 



.(3) 



.(3) 



(28) 



°1,2 u l,3 

The above can be inverted N times to give 

(3) (3) (3) (3) 

(3) C 2M C 0M-1 C 2,N-l C 0,N-2 



r (3)„(3) 
L 2,2 L 0,1 



'l.JV+l 



.(3) 

-1,JV-1 



J3) _ 
-l,N-2 



,(3) 
"1,1 



(3) (3) (3) (3) 

C 0,Af+l C 2,JV+2 c 0,AT+2 C 2,Af+3 



-(3) _ 
-l.JV+2 



J3) 

■"i.,N+Z 



(29) 



This equation with infinite continued fraction on the 
right-hand side can be solved numerically by minimizing 
the absolute value of the difference between the left- and 
right-hand sides. The equation has an infinite number of 
roots (corresponding to the QN spectrum), but for each 
TV the most stable root is different. In general, we have 
to use the N times inverted equation to find the iV-th 
QN mode. The requirement that the continued fraction 
be itself convergent allows us to limit its depth by some 
large value, always ensuring that an increase in this value 
does not change the final results within the desired pre- 
cision. To check the correctness of the obtained results 
we compared them with those obtained in with the 
help of WKB method 0. The WKB method is known 
to give very accurate results for higher values of I |25| . 

The complete numerical code for determining the QN 
spectrum was written in MATHEMATICA® , and is 
available from the last author upon request. 



IV. NUMERICAL RESULTS 

We now proceed to the presentation of our numerical 
results for the QN spectra of brane-localised fields with 
spin s = 0, 1/2 and 1 living on a higher-dimensional Kerr 
black-hole induced background. In what follows, we dis- 
cuss separately each species of fields and the dependence 
of its QN spectrum on the angular momentum of the 
black hole for arbitrary dimensionality of spacetime. 



A. Brane-Localised Scalar Fields 

Starting with the case of brane fields with spin s = 0, 
in Table I we display the fundamental modes (n — 0) of 
the QN spectrum for the partial mode (£ = 0, m = 0). 
Throughout this paper, unless explicitly stated, we will 
be focusing on the fundamental modes of the spectrum 
as they are found to dominate the signal, with the con- 
tribution of the second mode coming out to be about 1%, 
and that of the higher overtones even smaller. In Table I, 
the quasi-normal frequencies for the scalar field are pre- 
sented in terms of the angular momentum of the black 



Table I: QN frequencies loth for scalar fields, for the mode 
I — 0, m = and various values of the BH angular momentum 
a* = a/rh and bulk dimensionality D 



a, 


D = 5 


D = 6 


D — 7 


0.0 


0.27339-0.4109 li 


0.21040-0. 57167i 


0.07281-0.61990i 


0.5 


0. 24918-0. 35635i 


0.18233-0. 51895i 


0.04666-0.54481i 


1.0 


0. 19904-0. 26464i 


0.12312-0.42301i 


0.01619-0.40662i 


1.5 


0.15277-0.19716i 


0.06734-0.33503i 




2.0 


0. 11874-0. 15493i 


0.03337-0. 26202i 




2.5 


0.09556-0. 12834i 


0.01790-0.20889i 




3.0 


0.08032-0.11010i 


0.01081-0.17195i 




3.5 


0.06978-0.09604i 






4.0 


0.06167-0.08478i 






4.5 


0.05513-0.07574i 






5.0 


0.04978-0.06843i 






5.5 


0.04536-0.06241i 






6.0 


0.04166-0.05736i 







hole, measured in units of the black hole horizon radius 
r/i, and for the values D — 5, 6 and 7 of the dimensional- 
ity of spacetime - note that an upper bound exists on the 
value of the angular momentum parameter of the black 
hole given by a™ ax = (D — 2)/2 [26| and following from 
the assumption that the black hole was created by the 
collision of two particles with impact parameter b < r^; 
nevertheless, for completeness, here we extend our anal- 
ysis to higher values of a* thus allowing for black holes 
created by alternative mechanisms. From the entries of 
Table I, one easily observes that, as a* increases, both the 
real and imaginary part of the QN frequency decreases. 
The latter feature makes the damping time longer and 
the field oscillations on the brane longer-lived. On the 
other hand, for fixed a* and variable D, we find the same 
behaviour as in the case of a non-rotating black hole 
with the real part of the QN frequency being suppressed 
and the imaginary part enhanced, as D increases. 

As the reader notes, the study of the mode {I = 0, m = 
0) is incomplete; this is due to the very slow convergence 
of the Frobenius procedure for this particular mode as 
D and a become large. Although upon allowing suffi- 
cient computing time the missing quasi-normal modes 
could be found, the timescale for doing so does come out 
to be unrealistically large. For this reason, in Table II 
we present detailed numerical results for the next par- 
tial mode (£ = l,m = 0). These results confirm the 
monotonic suppression of both the real and imaginary 
part of the quasinormal frequency, as a* increases, thus 
pointing to the conclusion that this behaviour is inde- 
pendent of the particular multipole number. In addition, 
the more complete display of results available for this 
mode reveal some additional features of the spectrum. 
To start with, the increase in WR e found in [9j for the 
mode (£ = l,m = 0) as D goes from 5 to 6 is actu- 
ally followed by a decrease as D increases further, with 
this behaviour persisting for all values of a* ; on the other 
hand, the behaviour of wi m continues to be monotonically 
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Figure 1: Fundamental QN frequencies of a scalar field, as a 
function of a,, for the multipole mode 1=1, spacetime di- 
mensionality D — 5, and m = 1,0, —1 (depicted by triangles, 
rhombuses and squares, respectively). 



increasing as D takes larger values. 
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Figure 2: Quality factor for the scalar field multipole mode 
(t = 1, m = 0), and for D = 5, 6, ...9 (from top to bottom). 



One additional important feature emerging from the 
entries of Table II is the existence of asymptotic values 
for both WR e and uoim, as a» — ► oo. These asymptotic val- 
ues appear for all values of D, and for all scalar modes. 
In Fig. 2] we display the QN spectra for the scalar modes 
(I = 1, m = —1,0, 1) as a function of the black hole an- 
gular momentum parameter a*, for fixed dimensionality 
(D = 5). The presence of asymptotic values for both WR e 
and o;i m for all modes is obvious. Coming back to the en- 
tries of Table II, one may observe the potential existence 
of an asymptotic value for the imaginary part of the QN 
frequency also as D adopts large values: although LUi m 
clearly follows an increasing pattern, the step becomes 
gradually smaller as D increases - we will return to this 
point later on. 

Finally, a point to be addressed is the quality factor of 



the black hole as an oscillator, defined by 

The larger the above factor is, the better oscillator the 
corresponding field is. We find that the quality factor for 
scalar fields has a distinctly different behavior for differ- 
ent values of D: in the case of D = 5 there is a mono- 
tonic increase of the quality factor Q, as a function of 
the rotation parameter a, that renders the scalar fields 
better oscillators as the black hole rotates faster. How- 
ever, as D increases, the increase in Q is soon followed 
by a suppression, and, for large enough values of D, Q 
is monotonically decreasing. This behavior is depicted in 
Fig. |3 for the scalar mode (£ = l,m = 0). Although, 
in principle, the exact value of the quality factor is de- 
termined by all parameters of the black hole (Mbh, a ) 
and of the propagating field (s, i, m, n), it seems that, 
in the case of scalar fields and of the lower multipole 
modes, there are two main rival factors: the angular mo- 
mentum of the black hole, that tends to increase Q, and 
the spacetime dimensionality, that tends to decrease it. 
As is evident from Fig. |3 in all cases, the quality fac- 
tor reaches an asymptotic constant value as a* — > oo, in 
agreement with the discussion above. 

B. Brane-Localised Gauge Bosons 

We now turn to the case of gauge fields propagat- 
ing in the induced-on-the-brane background of a higher- 
dimensional rotating black hole. Our results are dis- 
played in Table III, and similarly to the case of scalar 
fields, they reveal the monotonic suppression of both the 
real and imaginary part of the QN frequency, as the an- 
gular momentum of the black hole increases. In addition, 
as a, oo, constant asymptotic values emerge for both 
WR e and wi m of the quasinormal frequency. This is obvi- 
ous from the results displayed in Table III for the mode 
(£= 1 , ttl — 0) , but it has been found to hold also for the 
modes m = ±1. 

When a* is kept fixed while D varies, a>R, e first in- 
creases as D goes from 5 to 6, in accordance with the 
generic behaviour found in 9], but as D increases fur- 
ther, a monotonic suppression takes over again. Also 
in agreement with the non-rotating case 0, the imagi- 
nary part of the QN frequency monotonically increases, 
as D grows, making the gauge field perturbations shorter- 
lived. The existence of a potential asymptotic value for 
LO\ m , as D — > oo, can again be inferred from the entries 
of Table III. 

The quality factor Q for gauge bosons is again strongly 
dependent on the value of the dimensionality of space- 
time D. Its behavior for the mode (£ = l,m — 0) and 
for the values D = 5,6, ...,9 is shown in Fig. As is 
obvious, for D = 5 we obtain again a sharp monotonic 
increase, as a* increases. For D = 6,7, we find, as in 
the case of scalar fields, a non-monotonic behavior with 
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Table II: QN frequencies tur h for scalar fields, for the mode i = 1, m = and various values of the BH angular momentum 
a* = a/r/t and bulk dimensionality D 







D 




5 




D = 


6 




D 


= 


7 




D = 


8 




D = 


9 


0.0 





75085- 




36387i 





81265-0 


50454i 





81678- 





61105i 





79233-0 


68002i 





76177-0 


71730i 


0.5 





68028- 





31425i 





74928-0 


45732i 





75414- 


o 


56410i 





72911-0 


63029i 





69952-0 


66372i 


1.0 





54711- 





23270i 





62853-0 


37640i 





63235- 





48047i 





60568-0 


53769i 





57925-0 


56198i 


1.5 


o 


43282- 


o 


17331i 


o 


52186-0 


31276i 





52101- 





40993i 


o 


49343-0 


45466i 





47184-0 


46989i 


2.0 





35010- 





13527i 





44102-0 


26836i 





43401- 





35623i 





40757-0 


38891i 





39096-0 


39791i 


2.5 





29112- 





11016i 





38023-0 


23672i 





36761- 





31440i 





34386-0 


33751i 





33138-0 


34279i 


3.0 





24803- 





09267i 





33336-0 


21315i 





31655- 





28076i 





29609-0 


29694i 





28670-0 


30011i 


3.5 





21554- 





07987i 





29619-0 


19478i 





27671- 





25307i 





25944-0 


26447i 





25230-0 


26642i 


4.0 





19033- 





07013i 





26599-0 


17991i 





24510- 





22993i 





23065-0 


23807i 





22512-0 


23930i 


4.5 





17025- 





06249i 





24096-0 


16751i 





21961- 





21035i 





20751-0 


21627i 





20317-0 


21707i 


5.0 





15392- 





05633i 





21988-0 


15690i 





19873- 





19364i 





18856-0 


19802i 





18508-0 


19854i 


5.5 





14041- 





05128i 





20190-0 


14766i 





18136- 





17925i 





17276-0 


18254i 





16994-0 


18289i 


6.0 





12904- 





04705i 





18640-0 


13950i 





16672- 





16675i 





15940-0 


16926i 





15708-0 


16951i 



Q 




Figure 3: Quality factor for the gauge field multipole mode 
(£ — 1, m = 0), and for D — 5, 6, ...9 (from top to bottom). 



an initial increase followed by a decrease. The deviation 
from the case of scalar fields arises for D > 7, where the 
monotonic suppression found in that case gives its place 
to a monotonic enhancement, although a marginal one. 



C. Brane-Localised Fermions 

We finally address the case of brane-localised spinor 
fields. In Tables IV and V, we present results for the 
two lowest multipole modes (£ — 1/2, m = 1/2) and 
(£ = 1/2, m = —1/2), respectively, for various values 
of the black hole angular momentum parameter and di- 
mensionality of spacetime. Unlike the previously studied 
cases of scalars and gauge fields, a deviation is observed 
in the behavior of these two modes, therefore we discuss 
both of them here separately. 

Starting from the mode (I = 1/2, m = 1/2), whose 
QN spectrum is presented in Table IV, we find that, in 
terms of a*, the generic behavior noticed for the other 
two species of fields also prevails here: both cjr and 



Wi m monotonically decrease as the angular momentum 
of the black hole increases, eventually reaching constant 
asymptotic values. On the other hand, for fixed a* and 
variable D, LUi m increases monotonically, as for the other 
two species of fields, with the only difference being the 
absence of an asymptotic value as D — > oo. A further 
deviation arises in the dependence of wr c on D with its 
exact behavior depending on the value of the black hole 
angular momentum: for low values of a* , wr c follows the 
traditional increase- followed-by-a-decrease pattern, while 
for high values of a* a monotonic decrease for all values 
of D is observed instead. 

Turning to the mode (£ = 1/2, m = —1/2) and to the 
entries of Table V, the deviation in the behavior of its 
QN spectrum is even more striking. As a* increases, its 
wr c and wi m monotonically decrease towards some con- 
stant asymptotic values, as expected, only for the values 
D = (5,6). As D increases further, the initial phase of 
decrease for WR e follows a phase of enhancement; toi m be- 
haves even more singularly by oscillating for D = 7 and 
monotonically increasing for D > 8. When a* is kept 
fixed and D varies, the previously observed behavior is 
followed only for low values of a*, with c<jr increasing 
until D = 6 and decreasing for greater values of D; for 
values of the angular momentum parameter higher than 
a* = 3, the enhancement phase is extended at least up 
to D = 7. Finally, wi m follows the traditional increas- 
ing pattern with D, however, as in the case of m — 1/2, 
there is no sign for the existence of an asymptotic value 
as D — > oo. 

The alert reader might have noticed that the deviations 
in the behavior of the QN spectrum for the aforemen- 
tioned spinor modes arise in the case of cither large a* 
or large D. In section IV. A. we commented on the very 
slow convergence of the Frobcnius procedure, for exactly 
these ranges of parameters, that failed to give us com- 
plete results for the lowest scalar mode (I = 0,m = 0). 
The same problem appeared also in the case of the spinor 
mode (£ = 1/2, m = —1/2), as it is obvious from Table 
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Table III: QN frequencies uorh for gauge fields, for the mode I = 1, m = and various values of the BH angular momentum 
a* = ajrh and bulk dimensionality D 



a* 




D 


= 5 




D 


= 


6 




D = 


7 




D 


= 


8 




D = 


9 


n n 


n 
u. 


O 1 UU 1 


fl ^1 74Qi 


n 

U 


c iS c !Q7 


n 
\j 


402 14i 


u. 






n 


559fi (S 


n 
•u. 


4797fli 


n 




,48451i 




n 
u. 


^9748 


fl 97^fi0i 

U.Z I OUOl 


U 


.J4UUO 


n 
•u- 


. OU J- 0£1 


u. 


^9789 D 


40874i 


n 

u 


.51222- 


n 

•u. 


4H1 RSi 


n 

u. 


c;nn^7 n 

UUUO ( U 


442731 


i n 

J. .u 








n 


*±ooo 1 


n 
•u. 


29121i 


n 
u. 


444Q7 fl 


QQfJCi4i 
OOUJ'll 


n 




n 
u. 




n 
\j 


4997fi fl 




1.5 


0. 


.34379- 


-0 14693i 


o 


.37737- 


•0. 


235741 


0. 


36855-0 


276191 





.35810- 


-0. 


9Q1 86i 

. £4 O J- QUI 


0. 


351 95-0 


99897i 




n 
u. 


970. 


nil ^07^ 


u 


^1 79^ 
.oil ZD- 


n 
-u. 


1 Q71 fii 

. iy / idi 


u. 


ouyuo-u. 


9^1 7*^i 
Zol i Ol 


n 
u 


^nnfi7 

.OUUD /- 


•u. 


9/1 ^7^ 


n 
u. 


9Q^fiQ n 
zyooy-u. 


9/1SH71 
.Z^ioU 1 1 


2.5 


0. 


.23272- 


-0.09115i 





.27183- 


•0. 


16975i 


0. 


.26399-0. 


19877i 





.25728- 


-0. 


20768i 


0. 


.25363-0. 


21086i 


3.0 


0. 


.19845- 


-0.07615i 





.23694- 


-0. 


14938i 


0. 


.22949-0. 


17365i 





.22410- 


•0. 


18042i 


0. 


22137-0 


18273i 


3.5 


0. 


17253- 


-0.06533i 





.20952- 


-0. 


13360i 


0. 


.20255-0. 


15396i 





.19817- 


•0. 


15919i 


0. 


.19609-0. 


16091i 


4.0 


0. 


.15238- 


-0.05718i 





.18751- 


■0 


12096i 


0. 


18105-0. 


13815i 





.17746- 


-0. 


14227i 


0. 


.17584-0. 


14358i 


4.5 


0. 


13633- 


-0.05083i 





.16950- 


■0 


11058i 


0. 


16356-0. 


12521i 





.16058- 


-0. 


1285H 


0. 


15931-0 


12952i 


5.0 


0. 


.12327- 


-0.04575i 





.15453- 


•0. 


10187i 


0. 


.14908-0. 


11443i 





.14659- 


-0. 


117121 





.14557-0. 


11792i 


5.5 


0. 


11245- 


-0.04159i 





.14190- 


•0. 


09445i 


0. 


.13693-0. 


10532i 


0. 


.13482- 


-0. 


10754i 


0. 


.13399-0. 


10818i 


6.0 


0. 


10335- 


-0.03812i 





.13113- 


•0. 


08804i 


0. 


.12658-0. 


09753i 





.12478- 


-0. 


09939i 


0. 


.12409-0. 


09991i 



Table IV: QN frequencies ur^ for spinor fields, for the mode I = 1/2, m = 1/2 and various values of the BH angular momentum 
a, = a/rn and bulk dimensionality D 



a* 




D 




5 




D 




6 




D = 


7 




D 




8 




D 




9 


0.0 


0. 


.44130- 


•0. 


35984i 





.45332- 


-0. 


51049i 


0. 


.41925-0. 


,65204i 





.33208- 


-0. 


78624i 





.31438- 


-1. 


28297i 


0.5 


0. 


.45308- 


•0. 


,29340i 





.45965- 


■0. 


43444i 


0. 


.41755-0 


,55905i 





.32135- 


-0. 


65497i 





.30694- 


-1 


11216i 


1.0 


0. 


.39498- 


•0. 


19854i 





.40166- 


■0. 


33780i 


0. 


.34902-0 


,44834i 





.25178- 


-0. 


49810i 





.23566- 


•0. 


95459i 


1.5 


0. 


.32503- 


•0. 


13271i 





.33293- 


-0. 


27005i 


0. 


.27090-0 


36141i 





.19483- 


-0. 


37575i 





.18243- 


•0. 


85971i 


2.0 


0. 


.26747- 


•0. 


09332i 





.27517- 


-0. 


22679i 


0. 


20951-0 


.2954H 





.15719- 


-0. 


29478i 





.14483- 


-0. 


,80778i 


2.5 


0. 


.22380- 


•0. 


06941i 





.22978- 


-0. 


1975H 


0. 


.16694-0. 


.2455H 





.13128- 


•0. 


2407H 





.11788- 


-0. 


77704i 


3.0 


0. 


.19081- 


•0. 


05414i 





.19420- 


•0. 


17589i 


0. 


13772-0 


,20815i 





.11251- 


-0. 


20284i 





.09839- 


•0. 


75720i 


3.5 


0. 


.16548- 


•0. 


04387i 





.16604- 


-0. 


15868i 


0. 


.11695-0. 


,17989i 





.09836- 


-0. 


17507i 





.08409- 


-0 


74360i 


4.0 


0. 


.14561- 


•0. 


03665i 





.14356- 


•0. 


14430i 


0. 


10155-0 


,15807i 





.08734- 


-0. 


15390i 





.07332- 


•0. 


,73388i 


4.5 


0. 


12969- 


•0. 


03137i 





.12551- 


-0. 


13196i 


0. 


.08972-0. 


,14082i 





.07852- 


-0. 


13726i 





.06498- 


-0 


72671i 


5.0 


0. 


.11671- 


•0. 


02740i 





.11091- 


-0. 


12125i 


0. 


.08036-0. 


,12689i 





.07131- 


-0. 


12386i 





.05835- 


-0 


72126i 


5.5 


0. 


.10594- 


•0. 


,02433i 





.09900- 


■0. 


11190i 


0. 


.07277-0. 


11543i 





.06531- 


•0. 


11283i 





.05296- 


•0. 


71701i 


6.0 


0. 


.09688- 


•0. 


,02192i 





.08918- 


•0. 


1037H 


0. 


.06649-0. 


,10585i 





.06024- 


-0. 


10360i 





.04849- 


•0. 


71363i 



V. Although numerical results were successfully derived 
for smaller values of either a* and/or D, it might be pos- 
sible that the looming bad convergence behavior of our 
numerical method contaminated these results, too. An- 
other reason for the observed modifications might be the 
fact that the equations for the QN modes corresponding 
to different Ps are entangled and, at times, it is extremely 
difficult to distinguish QN modes by their multipolc num- 
ber, especially in the case of ultra-spinning black holes. 
For the above reasons, the determined results for spinor 
fields for high values of either a» or D should be taken 
with caution. 

Finally, the quality factor for the spinor field mode 
(£ = 1/2, m = 1/2) closely follows the behavior observed 
in the case of scalar and gauge fields, as it is obvious from 
Fig. 0J For D = 5, Q has a monotonically increasing de- 
pendence on a* , that brings it to a value higher than for 
other species, thus rendering the fermions the best oscil- 
lators among the lowest multipole modes. For D > 6, the 
quality factor reaches a maximum value and then follows 
a decreasing pattern; for large values of D, the change in 



the value of Q, as a* varies, is only marginal. Turning 
to the mode (£ — 1/2, m = —1/2), we may easily see 
from Fig. that the quality factor follows a radically 
different behavior. Although the suppression with the 
dimensionality of spacetime is still present, the enhance- 
ment in terms of a* has disappeared. For D — (5,6), 
we can observe a monotonic decrease in Q, as a» grows, 
while for D = 7 the decrease in Q changes to an increase, 
after reaching a minimum at a* ss 2.5. The absence of 
robust results for this mode for D > 7 unfortunately 
does not allow us to check the dependence of its qual- 
ity factor on a* for larger values of D. Nevertheless, a 
clear conclusion can be safely drawn when the two lowest 
multipolc spinor modes are compared: the "co-rotating" 
modes with m = 1/2 are much better oscillators than the 
"contra-rotating" ones with m = —1/2. 
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Table V: QN frequencies loth for spinor fields, for the mode I = 1/2, m = — 1/2 and various values of the BH angular momentum 
a* = a/rh and bulk dimensionality D 





JJ = 


L) = 


L) — I 


L) = a 


L) — 9 


0.0 


0.44130-0.35984i 


0.45332-0.51049i 


0.41925-0.65204i 


0.33208-0.78624i 


0.31438-1. 28297i 


0.5 


0.35706-0.32705i 


0.37986-0.49522i 


0.34876-0.6624H 


0.25678-0.84279i 


0.17791-1.40707i 


1.0 


0.26937-0.25567i 


0.30417-0.43439i 


0.27267-0.63077i 


0.17717-0.87878i 


0.04019-1.48580i 


1.5 


0.20663-0. 198651 


0.25106-0.37925i 


0.21963-0.60756i 


0. 14808-0. 92662i 


- 


2.0 


0.16466-0.15934i 


0.21607-0.33850i 


0.18960-0.60173i 


0.16145-0.97336i 


- 


Z.o 


u.iooo f-u.ioiyoi 


u. iyz4o-u.ouyoyi 


u. i < oyo-u.ouoyii 






3.0 


0.11525-0.11217i 


0.17585-0.28853i 


0.18277-0.61502i 






3.5 


0.09988-0.09735i 


0. 16372-0. 27344i 


0.19531-0.62133i 






4.0 


0.08804-0.08588i 


0. 15462-0. 26239i 


0.21210-0.62451i 






4.5 


0.07865-0.07678i 


0. 14768-0. 25420i 


0.23044-0.62496i 






5.0 


0.07105-0.06939i 


0. 14236-0. 24804i 


0.24898-0.62352i 






5.5 


0.06477-0.06327i 


0. 13830-0. 24332i 








6.0 


0.05949-0.05814i 


0.13527-0. 23962i 









CONCLUSIONS 



1.5 



0.5. 




Figure 4: Quality factor for the spinor field multipole mode 
(I = 1/2, m = 1/2), and for D = 5,6, ...9 (from top to bot- 
tom). 



In this paper, we have performed a comprehensive 
study of the fundamental quasinormal modes of all Stan- 
dard Model fields propagating on a brane embedded in 
a higher-dimensional rotating black hole spacetime. The 
starting point of our analysis was the equations of mo- 
tion for fields with spin s = 0,1/2 and 1 propagating 
on the induced-on-the-brane background, and their re- 
writing to the form of a Schrodinger-like, wave equation. 
By using numerical techniques, based on the continued 
fraction method, the quasinormal spectra of all species of 
fields were derived as a function of the black hole angular 
momentum and dimensionality of spacetime. 

The dependence of the QN spectra on the angular mo- 
mentum of the black hole was found to be universal with 
both the real and (the absolute value of the) imaginary 
part of the QN frequency being suppressed as a* in- 
creases. As a result, brane-localised fields propagating on 
the induced 4-dimensional background are longer-lived 
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Figure 5: Quality factor for the spinor field multipole mode 
(I = 1/2, m = -1/2), and for D = 5, 6, 7 (from top to bot- 
tom). 
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Figure 6: Fundamental QN modes of fields with spin s 
0,1/2 and 1, for D = 6. 
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Figure 7: Fundamental QN modes of fields with spin s 
0, 1/2 and 1, for a* = 1. 



radically for modes with m < 0, whose quality factor 
is suppressed with D but also with a*, rendering them 
worse oscillators than the ones with m = or m > 0. 
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when the higher-dimensional black hole rotates faster. 
In addition, both WR e and wi m approach constant asymp- 
totic values, as a* — > oo. These results hold for all species 
of fields, with the only exception arising for spinor fields 
with m = —1/2 and for D > 7. Figure[S]shows the collec- 
tive behaviour of all SM fields QN spectra as a function 
of a*, for D = 6, and reveals the aforementioned spin- 
independent behaviour. 

In terms of the dimensionality of spacetime, our anal- 
ysis confirmed the increase in the imaginary part of the 
QN frequency of all SM fields, as D increases, observed 
in the case of a non-rotating black hole Q. The same 
behaviour holds for arbitrary angular momentum of the 
black hole, a result that points to the conclusion that an 
increase in the number of transverse-to-the-brane dimen- 
sions reduces the lifetime of the field perturbations on the 
brane, in both spherically- and axially-symmetric back- 
grounds. With the exception again of spinor fields, u)\ m 
tends again to a constant asymptotic value as D — > oo. 
This behavior is clearly shown in Fig. [7J where the QN 
spectra of all SM fields are shown for a* = 1. In the same 
figure, we also display the collective behavior of the real 
part of the QN frequency, cjr , as a function of D. This 
behavior is also found to be universal, with WR e increas- 
ing up to D = 6, in accordance with previously found 
results 9], but decreasing for higher values of D. 

The quality factor Q, finally, that determines the best 
oscillator among the different field perturbations, was 
found to depend both on properties of the field, like 
its spin s or multipole numbers (i,m), as well as on 
properties of the gravitational background, such as the 
dimensionality of spacetime and angular momentum of 
the black hole. The quality factor of modes with either 
m = or m > was found to decrease with D but 
(predominantly) increase with a*. As a result, for all of 
these modes, Q takes its largest possible value for D = 5 
and large a* - among the different species, the fermionic 
modes were found to assume the largest value of Q for 
these parameter values. The situation however changes 



Figure 8: A few higher overtones for the electromagnetic field 
mode [l = 3, m = 0), and D — 6. 

Being restricted by the standard model we did not 
compute quasi-normal modes of massive fields, yet, as 
was shown in |2?| for massive scalar and in |28| for mas- 
sive vector fields, the QN spectrum can show rather ex- 
otic behavior, for instance, infinitely long living modes. 
Also, throughout this paper, we have restricted our study 
to the fundamental overtones (n = 0) of the QN spectra 
due to their dominance in the signal. Nevertheless, higher 
overtones of all SM fields with spin s = 0, 1/2 and 1 were 
also looked at in the context of our analysis. Their QN 
spectra were derived for the range (0,77,.) of the angular 
momentum parameter a of the black hole with a step of 
r/,/8. In Fig. [SJ we display a few of the higher overtones 
for a brane-localised gauge field with (£ = 3, m = 0). For 
all higher overtones, it was found that an increase in a» 
leads again to the decrease of the absolute values of both 
the real and imaginary parts of their QN frequencies. 

We would like to finish our work with some impor- 
tant remarks. During our analysis, we made two sim- 
plifying assumptions. Firstly, our study was restricted 
to the case of higher-dimensional black holes with only 
one non-zero rotating parameter; that was due to the 
assumption that the colliding particles, that give rise to 
the black hole, were restricted to live on an infinitely-thin 
brane, and therefore had a non- vanishing impact param- 
eter only along an axis parallel to the brane. In a more 
realistic case, where the brane is assumed to have a finite 
thickness - generically of the order of the fundamental 
length - the colliding particles can have non-zero im- 
pact parameters along the bulk directions, and thus give 
rise to a black hole with more than one rotating param- 
eters. In more sophisticated brane-world models, that 
attack phenomenological problems like fast proton de- 
cay, large neutron-antineutron oscillations or FCNC, the 
quarks may be located at different 3-dimensional slices 
along the bulk and thus have, by construction, large sep- 
arations along the bulk directions [2{| . It seems therefore 
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imperative that a more comprehensive study of higher- 
dimensional black holes, with all possible components of 
the angular momentum taken into account, need to be 
done in the future (for scalar particles in D = 5, the gen- 
eral equation of motion has been derived and shown to be 
separable in pfflj L The second assumption made in our 
analysis was that the brane tension was much smaller 
than the mass of the produced black hole and thus it 
could be safely ignored. However, the "friction" between 
the black hole and the brane can alter the components 
of the angular momentum of the former |3lj |. and ex- 
cite brane degrees of freedom that could be interpreted 
as quasi-normal modes. The analysis of such a realistic 
set-up, that would study the change in the spectrum of 



the quasinormal frequencies found here due to the brane- 
black hole interaction, could also be the objective of a 
future work. 
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